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Abstract 

We use the N-quantum approach (NQA) to quantum field theory to construct 
a solution of the two flavor effective instanton induced 't Hooft interaction 
model valid for any temperature (T) and chemical potential (/i) beyond the 
mean field theory. The model contains only the qq channels. In constructing 
this solution we calculate the masses, widths and coupling constants of bound 
ex and resonant color 3 scalar diquark states. We find that the chemical po- 
tential induced by the interaction cancel exactly for all values of T in contrast 
to the Nambu-Jona-Lasinio model. Our method can also be extended to in- 
clude the qq channels in the model Lagrangian and we discuss how the NQA 
can be used to study the properties of bound diquark states at finite T and 
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I. INTRODUCTION 



In this letter we apply the N-quantum approach (NQA) to quantum field theory [[IJ to 
the effective instanton induced 't Hooft interaction Lagrangian || and construct a non- 
perturbative solution to the field theory beyond the mean field approximation for any 
temperature (T) and chemical potential (/x). Developed by Greenberg ||, NQA uses an 
expansion in on-shell asymptotic fields, known as the Haag expansion to solve the oper- 
ator equations of motion. The method is particularly well suited for bound state problems 
since it describes bound states with the same number of kinematical variables as in non- 
relativistic theories while maintaining explicit Lorentz invariance (for theories in free space). 
A summary of applications of NQA can be found in U . 

Recently we have extended the NQA to finite T and /x [H]. Our work was based on the 
thermal field dynamics formalism of Umezawa and collaborators || , and was motivated by an 
application of NQA to the BCS theory of superconductivity JTJ. Temperature and chemical 
potential were introduced simultaneously by applying a thermal Bogoliubov transformation 
to each of the fermionic or bosonic field operators in the Haag expansion. At the same time 
we defined a thermal vacuum state which is annihilated by each of the thermal annihilation 
operators. This takes care of the doubling of the Hilbert space which characterizes the 
formalism of thermal field dynamics. Furthermore, in the NQA one works only in the space 
of asymptotic field operators and therefore there is no need to specify the structure of this 
vacuum state. 

We have successfully tested our formalism on the Nambu-Jona-Lasinio (NJL) model 
H where the mean field (Hartree-Fock) results at finite T and /x are well established ||. 
These results were recovered using only the first term in the Haag expansion. In this work 
we take the second order terms in the Haag expansion into account to solve the 't Hooft 
model beyond the mean field theory. Although non-renormalizable and non-confining, the 
't Hooft interaction incorporates the symmetries of QCD and naturally provides a mechanism 
for chiral symmetry breaking and resolves the Ua(1) problem As is well known it also 



supports the formation of various kinds of condensates. Among them are the color singlet 
chiral condensate {(qq)) which breaks chiral symmetry dynamically and generates a mass to 
the fermions and the color anti-triplet scalar diquark condensate {(qq)) which triggers color 
superconductivity at finite density ]TT] . 



In this work we only include the qq channels in the model Lagrangian. Our solution is 
constructed by calculating the masses, widths and coupling constants of the a and resonant 
diquark states for any T and /x both in the Nambu-Goldstone and in the Wigner-Weyl 
phases of the model field theory. In the concluding section we discuss how our method can 
be used to study the properties of the bound diquark states in hot and dense matter. 

II. MODEL LAGRANGIAN AND THE NQA ANSATZ 

We start by defining the model Lagrangian which is the SU(2)f £g> SU(3) C version of the 
't Hooft interaction || with massless quarks and a chemical potential term. 

£ = z^(^-^ o7o )^ + g [(^*)V(^ 75 ^) 2 -(^7 5 *) 2 -(W*) 2 ] (2.1) 
= £njl - 9o [(^75*) 2 + (^m 2 } (2-2) 
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In the above Lagrangian ^fir^^, ^75 \l/ and \l/f\I/ terms correspond to the qq channels 
with quantum numbers of the a, tt, 7/ and mesons, respectively. The coupling constant 
go has the dimension of [mass] -2 . In Eq. ( |2.2| ) the 't Hooft Lagrangian is expressed as a sum 
of the NJL lagrangian and interactions in the 7/ and a channels. Note that /io is in general 
not the total chemical potential of the interacting many body system since the interaction 
may induce a contribution to the total chemical potential. 

For our purpose the equation of motion corresponding to Eq. ( |2.1| ) is most conveniently 
written in terms of colored quark fields in momentum space. For example, the charge 
conjugation symmetric form of the equation of motion for the up quark field may be written 
as 



H + ^lo)ap{u a )p{q) = -2g J d l p 1 d l p2d^p 3 5 4 (q -p x -p 2 - Pa) 



+ (75)a/?(75) 



(d b )e(pi),(d b ) p (p 2 ) {U a )p{pz)+ (4)e(pi),(4)/?(p2) {Ub)p{pz) 



(2-3) 



where [A, £?]_ = AB — BA. u and d are up and down quark fields with color and Dirac indices 
indicated in Roman and Greek letters, respectively. Summation convention for repeated 
indices is implied. 

We proceed to construct an ansatz for the 't Hooft model by Haag expanding each of 
the quark fields in Eq. ( |2.3|) to second order following |12|]. For example, the second order 
Haag expansion for the up quark field with color index 1 may be written as 



(«l)«(g) = K)a(g)+ / d 4 kd%5\q + k-b) \f a0 (k,b) : {u^ p{-k)a\b) : 



+d aP (k,b)(C- 1 )^ : \(d° 2 U-k)D° 2 (b) + (d*U-k)D° 5 (b) 



(2.4) 



where the symbol : : denotes normal ordering. The asymptotic fields in the Haag expansion 
are choosen to be in-coming fields which implies that the corresponding propagators obey 



the retarded boundary conditions. These fields are defined as u° a {k) 



u 



IN 



df{k)5{k 2 - m 2 ), a°(b) 



a 



IN 



(b)5(b 2 



and D° A (b) 



n iN 



(k)S(k 2 



m 



(b)5(b 2 



in 



D) 



In 



d° a (k) 

Eq. (|2.4j) o is a color singlet scalar field which transforms as while Da is a color anti- 
triplet scalar field transforming as ^C^j^Xa^ with A = 2, 5, 7 and C -1 = — C = —27270 
being the charge conjugation operator. 

In this work we only include the a and Da fields in the second order Haag expansion. 
Other composite fields may be included in the second order terms as shown, for example, 
in [|T2J. We emphasize that a and Da fields need not be interpreted as bound states. As we 
shall see below, in medium they can be either bound or resonant states. Thus, in order to 
avoid confusion with the bound diquark state |TT| we shall henceforth refer the color 3 qq 
state simply as the Da state. 

The asymptotic quark mass m will shortly be identified as the chiral gap while the a 
and Da masses m a and will be determined below as functions of T and /i. The Haag 
amplitudes f(k, b) and d(k, b) are interpreted as amplitudes for the creation of the a and 
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Da states, respectively. Corresponding expansion for (ui) a (q) can be found by using the 
relation (ui)(q) = u\(— (7)7°. Second order Haag expansions for the remaining quark fields 
in Eq. ( |2.3| ) can be constructed in a similar manner. 

We use the method developed in || to introduce temperature and chemical potential 
into Eq. ( [2.4| ). The on-shell fermionic field \1/ is extended to finite T and \x as follows.[] 

#°(g) = ty(q)5(q 2 -m 2 ) 

| TYIj 

= 7—7? — X) 6 (& s ) u (Qi - Wq) + d \~Q, s)v(-q, s)5(q° + U q 

{2ir) 6 uj q s=± l 

-> 7—7^ — X { a Q B (^ s ) ~ s Pq B \-Qi s ) u{q, s)5(q° - UJ q ) 



+ 



7,£> t (-g, s) - s5 q D(q, s) v(-q, s)5(q° + u q ) 



(2.5) 



Here B and are quasi-particle annihilation and quasi-hole creation operators, respec- 
tively, while D and are the annihilation operator for a quasi-anti-hole and creation 
operator for a quasi-anti-particle, respectively. These thermal operators satisfy the fermion 
anti-commutation relations and annihilate the interacting thermal vacuum state \Q(T,n)). 



B(q,s)\G(T,n)) — B(q,s)\Q(T, [m)) — D(q,s)\G(T, fi)) — D(q, s)\G(T, fM)) — 



(2.6) 



The T and \l dependent coefficients of the thermal Bogoliubov transformation are given by 

" 3 K±aO/(A; s T) 



a q = yjl — n~, (3 q = yjn~, 7 9 = ^Jl - n+ and 5 q = ^Jn+ where = [e^'/M + 1] 1 



are the Fermi distribution functions. 

The extension to finite T and \i for the on-shell bosonic field (fP proceeds through an 
analogous path. 

<f>°(q) = <P{q)5{q 2 -m 2 ) 
1 



2u. 



q 



a(q)5(q° - uo q ) + a\-q )6(q° + uo q 



2uo n 



jA(q) + d q A\q 



5(q° - u q ) 



+ 



d q A{-q) + CnA\-q) 5(q° + uj q ) 



(2.7) 



l + d 2 q 



The coefficients of the bosonic thermal Bogoliubov transformation are given by c q = 

and d q = \fnf where = [e^ - ^/^ 13 ^ — l]" 1 is the Bose distribution function. The thermal 
operators now obey the boson commutation relations and annihilate the thermal vaccum as 
in Eq. Q, 



Atf)\g(T,ii)) = Atf)\g(T,ri) = o 



(2.8) 



1 We use the conventions and normalizations of Itzykson and Zuber [13] 
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Thus our ansatz for the model Lagrangian Eq. fl2.1|) at finite T and fi is constructed 
in three steps. First, each of the colored quark fields in the equations of motion are Haag 
expanded to second order as in Eq. ( |2.4| ). Second, each of the asymptotic quark fields in 
the Haag expansion are subjected to a thermal Bogoliubov transformation of type shown 
in Eq. (|2.5| ). Finally, each of the asymptotic bosonic fields a and D\ are subjected to a 
thermal Bogoliubov transformation of type shown in Eq. (|2.7| ). The basic idea of NQA is to 
use the operator equation of motion Eq. ( |2.3| ) to solve for the Haag amplitudes f{k,b) and 
d(k, b) and thus construct a solution to the field theory defined by Eq. ( [2.1| ). 



III. THE GAP AND SELF CONSISTENCY EQUATIONS 

We now derive the gap equation and the self-consistency equations for the Haag am- 
plitudes using the operator equation of motion in the one-loop approximation. This ap- 
proximation corresponds to keeping only those terms involving a single contraction when 
renormal ordering the field operators with respect to the interacting vacuum \Q(T,/j,)). We 
begin by inserting the first order terms in the Haag expansion of the quark fields in Eq. ( ^.3| ) 
which is equivalent to invoking the mean field approximation. After renormal ordering and 
keeping only the linear terms in the field operators we obtain the following equation 

m = (24 + 4)g Q J d 3 p ^ (l - (3 2 p - 6%) (3.1) 

where d^p = d^p/(2ir) 3 . The factors of 24 and 4 represent contributions from the Hartree 
and Fock terms, respectively. In our approach both of these contributions arise naturally 
from the equation of motion. This gap equation determines the asymptotic quark mass m 
which is interpreted as the chiral gap. As in the NJL model when m(T, /i) = the model is 
in the Wigner-Weyl phase while if m(T, /i) ^ then it is in the Nambu-Goldstone phase. 

We find that the chemical potentials induced by the scalar and pseudoscalar Fock terms 
in the 't Hooft interaction cancel exactly.0 Therefore the total chemical potential, which 
appears in the Fermi and Bose distribution functions, is simply the bare chemical potential 
fio. This situation is in contrast with the NJL model where in the mean field approximation 
the Fock term induces a chemical potential proportional to the quark number density making 
it necessary to solve the gap equation self-consistently with the total chemical potential |H 



As a result, it becomes impossible to clearly define critical points in the T — plane |T4 
Eq ( |2.2|) shows that the 't Hooft Lagrangian can be written as a sum of the NJL Lagrangian 
and interactions in the a and the r( channels. Thus the addition of these channels leads to 
the cancellation of the chemical potential induced by the a and the ir channels. This is one 
important difference between the NJL and the 't Hooft model at finite \x in the mean field 
approximation. 

To go beyond the mean field theory we add the second order terms in the Haag expansion 
to our ansatz for the quark fields. Then, in addition to Eq. ( |3.1| ), we obtain the following self- 
consistency equations for the Haag amplitudes from the equation of motion in the one-loop 
approximation 



2 They are identified as coefficients of the 70 operator. 
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and 



(P-t-m)f(k,b) = 
2g J d 3 p^j(l-2/3. 



f(p~, b) + f(p + , b) + 75 /(p-, 6)75 + Kf(p + , b) l5 



+ 3Tr(/(p- 6) + /(p + , &)) + 3Tr( T5 /(p-, 6) + iJ(p + , 6)) 75 



/(p+ 6) + /(p", 6) + 75 /(p + , &)7s + 75/(p-, b) l5 



+ 3Tr(/(p+ 6) + f(p~, b)) + 3Tr( 75 /(p+, 6) + y 6 f(p-, b) 



75 



(p-l/k-m) d(k, b) 



m 
p—< 


y 


+ 


cr 


+ 




+ 




o ± = 





-ITjT/ 



d(p + ,b) +y 5 d(p + ,b)'j 5 



iT 



»-1TjX/ 



d(p , 6) + 7 5 c/(p , 6)75 



+ C' 11 d 1 (p~ ,b)C +i,C- 11 d 1 (p~,b)C 1 l5 



(3.2) 



(3.3) 



and 6 2 = m|, for Eq. (|3.3| ). These three dimensional linear integral equations for f(k,b) and 
d(k,b) can be solved exactly. 



IV. THE cr AMPLITUDE 



In the broken symmetry phase where m/0, the discrete C, P and T transformations 
constrain the structure of the a amplitude f(k, b) to be 



f G (k,b) 



1 



p — ft — m 



(A + BP) 



(-ft + m) 
2m 



(4.1) 



where A and B are real and constant diagonal matrices. The superscript NG indicates that 
we are in the Nambu-Goldstone phase of the model. The first factor in Eq. ( |4.1| ) is the 
quark propagator obeying the retarded boundary condition. At finite /1 it is necessary to 
choose a reference frame since Lorentz invariance is lost at finite density. We shall work in 
the rest frame of the bound states which is the appropriate frame to calculate their masses, 
widths and coupling constants as will be shown shortly. 

Upon inserting Eq. (|4. 1|) into the self-consistency equation for f(k, b) we find the follow- 
ing relation between A and B 



(A + Bp) (-ft + m) = (I f A + 2mB) (-ft + 



m) 



(4.2) 
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where 



If = 56go / d 4 p 5(p 2 — m 2 



2m 2 -b- P (i_o2_ r2 

(b - p + ie ) 2 - m 2 V P p > 



(4.3) 



with 6 M = (±m„, ). Following [12|] we trace Eq. ( |4.2| ) after multiplying by 1 and p to obtain 
two homogeneous equations for A and B. The solution for this set of equations exists when 
If = 1 which implies 



f NG (k,b) 



(A + 2mB) (-# + m) 
ft — ty — m 2m 



(4.4) 



It remains to determine the T and // dependent normalization constant A + 2mB. 

From the structure of the second order Haag expansion Eq. ( |2.4| ) we see that this nor- 
malization constant corresponds to the aqq coupling constant g aqq . This coupling constant, 
as well as the a mass and its width, can be determined by rewriting the condition If — 1 as 







2gc 



28 / d 4 p 5{p 2 - m 2 



2m2 ~ b -p h-^-S 2 

(b - p + ie ) 2 - m 2 V p 



(4.5) 



In the standard treatment of the NJL model the a mass is determined by the lowest pole 
of the a propagator in the random phase approximation (RPA) at zero three momentum. 
This pole is a solution of the equation 



= 1- 2GU a (m a ,0 



(4.6) 



where G is the coupling constant in the NJL Lagrangian and 11^ is the a polarization 
function |§. Since we are working in the rest frame of the a we can identify the term in 
square brackets in Eq. ( |4.5| ) as the a polarization function in the broken symmetry phase 



U^ G (b) = 28 J d 4 pd(p 2 - 



m 



2m 2 — b ■ p 



(l-/?-5 p 2 



(4.7) 



(b — p + ie ) 2 — m 2 

For 6 M = (b , ) and 6q > 4m 2 , the real and imaginary parts of Eq. (|4.7|) are given byf\ 



7 r°° 
Re r£ G (6 ) = — P / dx 



Im n^ G (6 



7T- 

7 1 

7T&o 
rNG 



{x 2 -m 2 fl 2 \ 



- 

i 3/2 



1 7,2 
4°0 



4 



b — m z 



i -i 



,(x+lS)/(k B T) 



e (M +M )/ (fcsT ) 



(4.8) 
(4.9) 



The real part of 11^ is equivalent to the one obtained by Asakawa and Yazaki in Eq. (30) 



of [|1J] except for the factor in front of the principal value integral. They obtain a factor of 



The real part of n^ G is most easily determined by calculating the imaginary part and relating 



the result to the real part as shown in 15]. 
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Y whereas our factor is 7.[] Note that Re n^ G is a function of 6q and is independent of the 
choice of the b ± = (±6 ,0) frames. 

The a mass can be obtained directly from Eq. ( |4.5| ) by noting that the gap equation 
Eq. can be rewritten as 



1 = 28^o / drp 5(p z - m 



2 m 2 a -2b-p 



(4.10) 



Using this expression in Eq. 
width of the a state is given by 



(b — p) 2 — m 2 

implies m a = 2m for all T and fi. Moreover, since the 



r„ 



i 



Urn EL 



(4.11) 



it follows immediately from Eq. ([4.9|) that the a has vanishing width in the broken symmetry 
phase. 

Finally, in the RPA the aqq coupling constant is given by |I 



2 

9 aqq 



dRe IT 

db 2 



(4.12) 



fi 2 =mj 



which fixes the normalization of the a amplitude in the broken symmetry phase. Thus, we 
have solved for the Haag amplitude of the a state which can be written as 



f G (k,b) 



g*qq 



-I/t + m) 



p — ft — m 2m 



(4.13) 



and obtained its mass, width and coupling constant. In the broken symmetry phase the a is 
a stable bound state with zero width and a mass of exactly twice the chiral gap for all values 
of T and /i. This conclusion is in agreement with the one obtained by Zhuang, Hiifner and 
Klevansky in the NJL model |16[ . 

To obtain the a amplitude in the chirally symmetric phase we repeat the same exercise 
elucidated above but with m = 0. The structure of the Haag amplitude now becomes 



f BW (k,b) 



g*qq 



(4.14) 



while the real and imaginary parts of the extracted a polarization function are found to be 



Ren^ w (6) = ^P 



7T- 



dx 



x 



x 2 — \m 2 a 
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87T 



e {\m a -P F )/{k B T) 



,{x-P F )/(k B T) 



e {\m a +P F )/{k B T) 



{x+P F )/{k B T) 



-1 



(4.15) 
(4.16) 



4 Note that their coupling constant G is normalized to G = 4 ^ t " 1 g where g is the coupling constant 
appearing in the NJL Lagrangian. 



S 



where Pp is the Fermi momentum. In this case the a mass must be determined by solving 
the equation 



= l-2g ReU^(b) 



(4.17) 



for each T and Pp. This mass determines the width and coupling constant through 
Eqs. OCT]) and (gig) . 

The a state in the Wigner-Weyl phase is in general not a stable bound state. For T/0 
it is a resonant state with a finite width which can decay into a quark-anti-quark pair. 
However, when T = and Pp ^ the decay width vanishes for m a < 2Pp and the a 
becomes a stable bound state. This is a consequence of Pauli blocking which allows the a 
to decay only when m a > 2P F [IH . 



V. THE D A AMPLITUDE 



Determination of the Da amplitude, mass, widths and coupling constant is analogous 
to that of the a case. In the Nambu-Goldstone phase the Da amplitude in the broken 
symmetry phase is given by 



d NG (k,b) 



NG 
iJDqq 



+ m) 



P — $ — m '° 2m 
where gp, qq is the Da coupling constant. The extracted polarization function is 



-b ■ p 



. dp5(p -m )— — . 
7r J [o — p + le'Y 



l-2%9(p )- 25 2 p e(-p ) 



(5.1) 



(5.2) 



Similarly, the Da amplitude and the polarization function in the Wigner-Weyl phase are 

,BW 



d BW (k,b) 



and 



n| w (fc) = ^ I d*p s(p 



7T° 



1 



—b ■ p 


(b — p + ie°y 

,{p-li)/(k B T) , 



9pqq 



(5.3) 



0(po) 



,(p+A»)/(fcflT) 



+ 1 



0(-p o ; 



(5.4) 



Note that the Da polarization functions Eq. (|5.2j ) and (|5.4|) are not symmetric under the 
interchange of particle and anti-particle Fermi distribution functions due to the presence of 
the Heaviside functions 0{p o ) and 9(— p ). As a consequence there are two possible branches 
of the Da state. For example the real part of II^ G on the positive mass shell b + = (+J71d, ) 
is 



ReU N D G (b 



2 I r°° 

-P/ <£r , . 

7T Z [ Jm X z — jm D 



in- _ f°° , x\J x 1 — m? 



1^2 



+ 2P 



2 r 



rap, -\-2x 



rri£) — 2x 



,(x-n)/(k B T) 



,(x+n)/(k B T) + 



(5.5) 



while on the negative mass shell b = (— m_o,0) the corresponding result is 



RenS G (n = A p| t '* J -2/ <fa 





e (x-n)/(k B T) + 1 



1 -1 



(5.6) 



The origin of the two branches is due to the unequal treatment of particles and anti-particles 
at finite density. When \i = 0, the particles and anti-particles are described by the same 
Fermi distribution function and thus the Da polarization functions become invariant under 
the b + and b~ interchange leading to a single branch. 

VI. a AND D A MASSES, WIDTHS AND COUPLING CONSTANTS 

We present results for the a and Da masses, widths and coupling constants for zero 
temperature and finite density. Results for other combinations of T and fi will be given 
elsewhere together with relevant technical details of this work. We use the non-covariant 
regularization and simply cut off the divergent integrals using a cut off parameter A. The 
other input parameter is the coupling constant go- We arbitrary choose A = 0.61 GeV 
and go = 5.01 GeV~ 2 which yields a dynamical quark mass of m = 0.305 GeV « \M N at 
(T, /i) = (0, 0), where M N is the nucleon mass.Q 

In Figure la we plot the a mass as a function of Fermi momentum. In the broken 
symmetry phase of the model m CT = 2m where m is dynamical quark mass. We see that 
ma decreases continuously with increasing Pp and vanishes at the critical Fermi momentum 
(-Pf)c °f about 0.3 GeV. The phase transition is second order like in the NJL model. Beyond 
(-Pp)c the a mass starts to increase monotonically and by Pp = 0.44 GeV its value is equal 
to the value in free space. For all values of Pp we find that m a < 2Pp and therefore it has 
vanishing width. Hence we find the a to be a stable bound state at finite density and zero 
temperature in this model. 

The behaviour of the density dependent a mass is reflected in the aqq coupling constant 
shown in Figure lb as a function of Pp. When Pp = the value of g^ qq /4iT is about 
0.342 and it continuously increases with density until the chiral symmetry restoration point 
where it assumes its maximum value of 0.636. Thus as the a mass decreases its coupling 
constant increases indicating the strengthening of the qq binding. Above (Pf) c , the coupling 
constant starts to decrease with the increasing a mass. However, when the a mass becomes 
comparable to its value in free space the coupling constant is about 50% larger than the free 
space value. Hence the a mass is much more sensitive to changes in its coupling constant in 
the Wigner-Weyl phase than in the Nambu-Goldstone phase. 

We find that there are no solutions to the equation determining the D A mass in the 
b~ frame with our choice of input parameters. The reason is the rather small value of the 
regulator A. We can find solutions for the Da mass in this branch if we cut off the divergent 



5 These are the same parameters used in ||. 
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integrals in the Da polarization function independently of A. However, we choose to work 
with only one regulator as in previous works of the NJL model at finite T and \i ||. With 
a single regulator we find solutions for the diquark mass in the b + frame but only for Pp > 
0.06 GeV. Thus there are no Da states in free space and for T = its formation occurs only 
above a critical density the value of which depends on the input parameters. This implies 
that below this critical density there is no contribution to the solution of the 't Hooft model 
from the Da term in the Haag expansion. We also find solutions for the Da width indicating 
that this is a resonant state in contrast to the a state. The Da mass and its width are shown 
in Figure 2a. 

In the broken symmetry mode the Da mass peaks around Pp = 0.2 GeV and decreases 
monotonically until the phase transition point. Beyond this point it increases almost linearly 
as in the a case. We find that nip, > m a for all values of density in which there is a resonant 
Da state. However, the Da coupling constant is one order of magnitude smaller than that 
of the a state as shown in Figure 2b. In fact, gp> qq does not become appreciable until Pp 
~ 0.2 GeV and continuously increases with increasing density. Since gp> qq normalizes the 
Da amplitude this means that the contribution from the Da term in the Haag expansion 
becomes important with increasing density but it will always be a correction to the dominant 
a term. 



VII. CONCLUSION 

In this work we have constructed a solution to the 't Hooft interaction Lagrangian valid 
for finite temperature and chemical potential beyond the mean field approximation. To our 
knowledge this is the first time that such a solution has been presented in the literature. In 
order to construct this solution we have used the NQA to quantum field theory involving 
the second order Haag expansion in the one-loop approximation. We have only included the 
color 1 a and color 3 Da composite fields in our ansatz to the 't Hooft model. An improved 
solution may be obtained by adding all possible channels, both color singlet and non-singlet, 
allowed by the model Lagrangian. 

Our solution was obtained by solving the self-consistency equations for the a and Da 
Haag amplitudes. This has been accomplished by calculating the masses, widths and cou- 
pling constants of the a and the Da states. At finite density and zero temperature we 
find that the a is a stable bound state with a vanishing width while the Da is a resonant 
state with a finite width. The latter state is resonant because we have only included the 
qq channels in our model Lagrangian. Thus our model provides no mechanism to create a 
bound Da state. These properties of the a and the Da states can also be calculated using 
the standard Hartree-Fock-RPA approach. We stress that the advantage of using the NQA 
over the Hartree-Fock-RPA approach is that one can explicitly obtain a solution to the 
model field theory. 

It would be interesting to apply our method to study the properties of a bound diquark 



state [|TTJ at finite T and \x. This may be done by adding a color 3 qq channel to our model 
and repeat the calculation presented here using the same ansatz. A vanishing diquark width 
would indicate a bound state just as in the case of the a shown in Figure la. In addition we 
can determine the color superconducting gap by calculating the off-diagonal Hamiltonian 
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and demanding that it vanishes. We hope to report on the outcome of this investigation in 
the near future. 
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FIGURE CAPTIONS 



FIGURE 1. a mass and aqq coupling constant as functions of Fermi momentum P F 
for T = 0. Input parameters are A = 0.61 GeV and g = 5.01 GeV~ 2 . The critical Fermi 
momentum is (-Pf)c = 0.301 GeV. (a) m a vs. P F . The a is a stable bound state and has 
vanishing width, (b) g^ qq /A7T. 

FIGURE 2. Color 3 scalar diquark mass and Dqq coupling constant evaluated in the 
frame = (+m£>,0 ) as functions of Fermi momentum P F for T = 0. Input parameters 
are the same as in Figure I. (a) m D vs Pp. The diquark state does not exists for P F < 
0.06 GeV. It is a resonant state and its width T is also shown in the figure bounded by the 
dashed line, (b) g 2 Dqq /A-n. 
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